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I.  INTRODUCTION 

This  paper  compares  four  major  schemes  used  for  forecasting 
demand  to  be  used  as  -input  to  an  inventory  model  sc  that  "optimum" 
stockage  levels  can  be  obtained.  Although  the  inventory  model  used 
(the  classical  periodic  review,  order  up  to  S  ,  infinite  horizon 
model  with  complete  backordering)  assumes  a  known  stationary  demand 
distribution,  it  is  seldom  true  that  demand  is  in  reality  stationary. 
This  study  assumes  demand  is  known  to  be  Poisson,  but  the  mean  de¬ 
mand  in  each  period  is  unknown.  Maximum  likelihood  (sample  mean), 
exponential  smoothing,  standard  Bayes  and  adaptive  Bayes  forecasting 
schemes  are  compared,  via  Monte  Carlo  simulation  in  situations  when 
the  mean  (1)  is' actually  stationary,  (2)  has  a  long  term  trend  and 
(3)  incurs  "shock"  changes  at  random  times.  The  criterion  for  com¬ 
parison  is  the  sum  of  the  total  discounted  costs. 

II .  INVENTORY  MODEL 

Denoting  the  holding  penalty  by  h  ($  per  unit  on  shelf  at 
the  end  of  a  period),  the  backorder  penalty  by  ir  ($  per  unit  in 
backorder  at  the  end  of  a  period,  the  unit  put  chase  cost  by  c 
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($  par  unit),  the  cumulative  distribution  function  of  the  demand 
par  period  by  F(*)  and  assuming  no  fixed  cost  of  ordering  and 
zero  leadtime,  the  policy  which  minimizes  the  total  sum  of  dis¬ 
counted  costs  over  an  infinite  horizon  with  discounting  factor  a 
is  given  by  (see,  for  example,  [5],  p.  386  or  [6]), 

"At  the  beginning  of  each  period  order 
max  [S*-  -  z ,  0]  , 

where 

z  ■  on  hand  inventory  before  ordering 
and 

(1)  F(S*)  -  [ir  -  c(l-a)  ]/ (h+ir) 

If  F(* )  is  discrete,  the  equality  in  (1)  becomes  a  "just  greater 
or  equal  to"  inequality,  that  is,  S*  is  the  smallest  integer  value 
for  which  F(S*)  is  Just  greater  or  equal  to  the  right  hand  side. 


Further,  it  can  also  be  shown  that  for  a  finite  horizon  of 

n  periods,  the  optimal  policy  is  to  order  up  to  a  single  critical 

value,  S  ,  where 
n 


S1  -  S2  - 


<  S  .  <  S  <  S* 
n-1  n 


and  S*  satisfies  Equation  (1).  Practice  has  shown  that  convergence 
to  S*  is  quite  rapid,  especially  when  the  cost  of  a  unit,  c  ,  is 
relatively  small. 

Equation  (1)  is  the  expression  used  in  the  simulator  to  de¬ 
termine  the  best  inventory  level  S*  for  each  period. 


III.  FORECASTING  SCHEMES  TESTED 

It  is  assumed  that  the  demand  per  period  is  Poisson,  with 
mean  A  ,  but  that  X  is  unknown,  and  must  be  estimated  from  past 
data.  Four  major  procedures  (actually  five  procedures  In  all)  are 
considered  in  this  study:  maximum  likelihood,  exponential  smoothing 
(and  exponential  smoothing  with  trend),  standard  Bayes  and  adaptive 
Bayes.  These  are  discussed  m  more  detail  below. 
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Maximum  Likelihood  (ML) 

It  ia  easy  to  show  that  the  maximum  likelihood  estimator  for 
the  mean  of  a  Poisson  random  variable  ia  the  arithmetic  mean  of  all 
past  observations  (see,  for  example,  [2],  pp.  12-13).  Thus,  this 
scheme  is  extremely  easy  to  use  and  furthermore  would  be  the  most 
familiar  to  those  with  only  a  passing  acquaintance  of  statistics. 

Exponential  Smoothing  (ES) 

The  simplest  of  the  exponential  smoothing  routines  is  given 
as  (see  [2],  p.  14,  [1],  p.  101  or  [4],  p.  B-128) 

(2)  Xn  -  axfl_1  +  (l-a)*^ 

where 

A 

1  ■  the  estimated  mean  for  period  n 

a  •»  smoothing  constant  (Olail)1 

x  ■  the  demand  observed  in  period  n  . 
n 

In  order  to  start  the  procedure,  we  use  the  Initial  condition 

A 

A2  -  (x2+x1)/2  , 
thus  (2)  ia  valid  for  n  >  3  . 


Exponential  Smoothing  With  Trend  (ESI) 


If  a  long  term  trend  is  expected,  better  results  are  gen¬ 
erally  obtained  by  including  a  trend  estimator  explicitly  in  the 
forecasting  scheme.  The  forecast,  *n  ,  now  becomes  (see  [2],  pp. 
14-15  or  [4],  p.  B-128) 


(3) 


X 

n 


n  n 


> 


used. 


1 


For  all  runs  throughout  this  study,  a  value1  of  0.3  was 


-  3  - 


where 
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* 

n 


b 


ax  +  (1-a)  .■#  ) 

n  n-l  n-l 

b(£  . )  +  (l-b)3?  . 

n  n-l  1  n-l 

2 

trend  smoothing  constant  (0<b<l) 


and  the  Initial  conditions  used  are 


i 


making  (3)  valid  for  n  >  2  . 


Standard  Bayes  (£B) 

For  the  Bayes  schemes  it  is  necessary  .o  assume  a  prior, 
distribution  on  X  ,  since  under  Bayesian  analysis,  A  is  not 
assumed  to  be  a  fixed  unknown  constant  but  a  random  variable,  The 
prior  distribution  assumed  for  A  is  an  exponential  distribution 

A 

with  a  known  mean  9  .  The  estimate  A^  is,  for  this  standard 
Bayes  analysis,  the  expected  value  of  the  posterior  distribution  , 
for  A  given  the  sum  of  the  past  n-l  observations  which  we 
denote  by  tR_^  •  The  posterior  distribution  given  t  ^  turns 
out  to  be  gamma  with  parameters  n  -  1  +  1/9  and  cn_^  +  1  (Bee 
[2],  pp.  9-10)  so  that  the  expected  value  and  hence  A^  is  given  by 

(4)  An  -  (t  +l)/(n-l+l/6)  . 

n  n-i 

h'ote  that  for  n  ■  1  and  =  0  ,  A^  is  just  the  mean  of  the  prior, 

namely  0  ,  so  that  (4)  holds  for  n  >  1  .  Also,  the  limit  of  (4) 
as  n  »  is  the  arithmetic  mean  so  that  for  large  n  ,  A^  giveh  by 
(4)  is  approximately  equal  to  the  maximum  likelihood  forecast. 


2 


] 


For  all  runs  throughout  this  study,  a  value  of  0.3  was  used. 
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Adaptive  Bayes  (AB) 


The  adaptive  Bayes  procedure  (see,  for  example,  Zachs  [7]) 
instead  of  using  the  expected  value  of  the  posterior  distribution 
as  the  Poipson  A'  in  the  inventory  equation  (Equation  (1)),  actu¬ 
ally  utilizes  the  marginal  distribution  of  demand  for  the  nth 
period  given  the  sum  of  the  past  n  -  1  periods  of  demand.  This 
distribution  function  is  used  for  F(0  in  (1)  and  S*  determined 

accordingly.  Denoting  this  marginal  distribution  as  g  (■)  and 
■  n 

its  cumulative  distribution  by  Gn(*)»  then  (1)  becomes 


(5) 


G  (S*)  >  [tt  -  c(l-cO]/(h+iT)  . 


It  turns  out  (see  [2],  pp.  10-12  or  [7])  that  G  (•)  is  a  negative 
binomial  with  parameters 


t1  .  +  1  and  tl/(n+l/6)]  90  that 

n-x 


S* 


S*  ■ 

(6)  G  (S*)  -  £  g  (x| t  f)  -  £ 

"  ,  x-0  n  x-0 


/x+tn-l]  /n -1+1 /el tn_1  V  1  V 

^  x  y^n+l/ey  ^  rH-l/ey 


IV.  SIMULATOR 

! 

Poisson  demand  is  generated  via  a  Monte  Carlo  procedure  with 
the  mear.  A  being  an  input  parameter.  In  addition,  a  value  of  G 
is  required  as  input  in  order  to  operate  the  Bayes  schemes.  The 
user  has, the  option  of  requiring  A  to  remain  constant  (stationary 

I 

case),  having  A  increase  by  any  specified  amount  each  period  (long 

1 

term  trend  case)  cr  having  A  "jump"  by  a  specified  amount  at  spec¬ 
ified  times  (shock  case).  The  trend,  and  jump  values  are  also  input 
parameters.' 

The  cost  parameters  (c,  h,  tt)  as  well  as  the  discount  factor 
a  and  initial  on  hand  inventory  must  also  be  inputted.  The  user 
also  specifies  the  number  of  periods  to  bs  run  and  the  number  of 
times  a  given  run  is  to  be  replicated  (replications  refer  to  keeping 
all  input  parameters  fixed,  but  generating  different  sets  of  random 
deviates  for  the  demand  stream  thus  allowing  for  statistical  analyses 
on  the  output).  An  initialization  period  must  also  be  specified 
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during  which  forecasts  are  made  and  stockage  levels  calculated  but 
resulting  costs  are  not  Included  in  the  discounted  cost  streams, 
thereby  removing  start-up  effects  from  the  cost  comparisons.  The 
program  also  calls  for  the  specifications  of  a  pre-initialization 
period  in  which  demand  values  are  generated  prior  to  calculating 
the  first  S*  so  that  every  forecasting  scheme  will  have  the  same 
number  of  observations  prior  to  making  its  first  forecast. 

The  user  must  specify  which  forecasting  schemes  are  to  be 
employed.  All  or  any  combinations  of  the  five  schemes  discussed 
in  Section  III  can  be  utilized.  When  more  than  one  scheme  is  called 
for  on  a  run  set  (a  run  set  refers  to  the  set  of  replications  for 
a  fixed  set  of  input  parameters)  the  same  demand  stream  is  used  for 
each  scheme.  That  is,  demand  is  generated  via  Monte  Carlo,  the 
first  scheme  is  used  for  forecasting  and  stockage  calculations  and 
costs  are  computed.  Then,  on  the  same  set  of  demand,  the  next  scheme 
is  used,  etc.  Thus,  one  replication  is  generated  for  each  scheme. 
After  the  last  scheme  Is  used,  a  new  demand  stream  is  generated  and 
the  process  repeated  to  generate  another  replication  for  each  scheme. 
The  procedure  continues  until  the  desired  number  of  replications  are 
obtained.  By  using  the  same  demand  stream  for  each  forecasting 
scheme  on  each  replication,  the  cost  differences  are  due  solely  to 
the  effect  of  forecasting  and  hence  the  variation  among  the  average 
discounted  costs  for  the  various  forecasting  schemes  is  kept  as  low 
as  possible.  Also,  the  cost  results  lend  themselves  nicely  to  sta¬ 
tistical  analysis  by  a  paired-sample  t  test. 

The  simulator  puts  out  various  quantities  of  interest,  the 
major  being  the  average  discounted  expected  costs  per  period,  that 
is,  the  total  sum  of  discounted  expected  costs  for  the  number  of 
periods  run  divided  by  the  number  of  periods.  Also,  average  on  shelf 
inventory,  average  number  in  backorder  and  the  demand  stream  generated 
as  well  as  the  resulting  S*'s  calculated  for  each  forecasting  scheme 
arc  available.  The  simulator  is  written  in  FORTRAN  IV  and  had  been 
run  on  the  IBM/360/50. 
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V.  STATISTICAL  TESTS 

A  chi-square  goodness  of  fit  test  (see,  for  example,  [3], 

pp.  285-6)  was  performed  on  the  Poisson  demand  generator.  A  sample 

of  180  periods  demand  was  generated  with  A  *=  5  .  Nine  class  inter- 

2 

vals  were  used  and  the  calculated  x  test  statistic  was  8.6448. 

2 

Prom  x  tables  (see,  for  example,  [3],  p.  318),  a  5%  critical  value 
for  8  degrees  of  freedom  is  15.507  and  a  10%  critical  value  is  13.362. 
Thus,  there  is  no  evidence  to  reject  the  hypothesis  that  the  generated 
data  came  from  a  Poisson  distribution  with  A  *  5  . 

To  determine  if  significant  differences  can  be  found  in  ex¬ 
pected  discounted  costs  when  using  different  forecasting  schemes,  a 
paired-sample  t  test  was  employed  (see,  for  example,  [3],  p.  269). 
This  test  utilizes  che  paired  differences  in  the  costs  between  any 
two  schemes  and  tests  if  the  average  difference  (over  all  replications 
of  a  given  run  set)  is  significantly  different  from  zero. 

VI .  RESULTS 

Prior  to  investigating  the  effect  of  forecasting  schemes  on 
inventory  costs,  two  preliminary  tests  were  made  to  determine  (1)  the 
sensitivity  of  the  assumption  on  the  prior  mean  (6)  in  the  Bayes 
schemes  and  (2)  the  effect  of  the  discount  factor  (a).  Next,  three 
major  cases  were  investigated  to  study  the  effect  of  forecasting 
schemes  on  inventory  costs  under  (1)  stationary  demand,  (2)  long  term 
trend  and  (3)  shock  changes  in  mean  demand.  Finally,  an  analysis  of 
maximum  likelihood  versus  the  Bayes  schemes  for  low  mean  demand  was 
performed . 

In  performing  the  simulations,  decisions  must  be  made  as  to 
the  set  of  costs  used.  These  cost  parameters  are  related  to  each 
other  in  the  following  way.  The  right  hand  side  of  Equation  (1)  (or 
(5))  is  a  particular  fraction  depending  on  the  values  of  it,  h  and 
c  ,  assuming  Ci  is  fixed.  Let  this  fraction  be  denoted  by  f  , 
that  is, 

(7)  f  -  U  -  c(l-a)]/(h+Tr)  . 
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Tables  2,  3  and  4  give  the  results  of  the  test  for  sensitivity 

4 

of  the  Bayes  scheme  to  8  .  Twenty  replications  of  runs  of  40  periods 
each  were  performed  for  X  stationary  at  5,  f  -  0.8  and  a  •>  3/4  . 

Table  2  shows  the  mean  and  standard  deviations  of  the  average  dis¬ 
counted  costs.  Table  3  gives  the  percent  difference  between  6  ■  5 
and  1  and  8  ■  5  and  9  for  each  Bayes  scheme,  when  the  true  mean  is 
5  and  the  demand  process  is  stationary.  The  largest  percent  difference 
in  Table  3  is  less  than  1.5%.  Table  4  shows  that  when  a  paired  t  test 
is  run  on  the  differences,  no  statistical  significance  at  either  the 
5%  or  1%  level  is  noted.  Thus,  for  all  other  runs  (except  the  final 
run  sets)  6  was  always  set  equal  to  X  and  it  was  felt  that  this 
would  not  provide  a  significant  advantage  to  the  Bayes  schemes  over 
Che  non-Bayes  schemes. 

Tables  5  and  6  give  the  results  of  run  set  2  to  determine  the 
sensitivity  to  a  assumptions.  The  cost  ratio  f  was  held  at  0.8, 
while  a  and  X  were  set  at  3/4  and  1,  and  5  and  1  +  0.5  respec¬ 
tively  (1+0.5  indicates  X  starts  at  1  and  increases  by  0.5  per 
period).  Again,  20  replications  of  runs  of  40  periods  were  performed. 
From  Table  5  we  see  that  the  relative  order  of  average  cost  per  period 
is  preserved  across  schemes  for  the  different  a's  .  Also,  from  Table 
6,  the  percent  differences  are  generally  the  same  order  of  magnitude. 

It  is  noted,  however,  that  when  discounting  takes  place  (a^l)  a  large 
reduction  in  the  average  cost  per  period  results.  This  is  due  to  the 
fact  that  the  cost  per  period  is  essentially  discounted  to  zero  for 
the  later  periods.  These  later  periods  with  zero  or  nearly  zero  period 
costs  are  averaged  in  with  the  earlier  non-zero  costs,  resulting  in  a 
greatly  reduced  average  period  cost.  Since  this  is  effectively  throwing 
out  (or  at  least  greatly  diminishing  the  influence  of)  later  period 
costs,  it  was  concluded  that  a  discount  factor  of  one  should  be  used 
which  would  allow  all  periods  to  be  weighted  equally.  (By  using  an 
a  ■  1  ,  we  are  actually  comparing  the  expected  costs  per  period  (cost 


^These  40  periods  are  over  and  above  10  periods  used  for 
initialization  plus  5  periods  used  for  pre-initialization. 
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TABLE  2 


RUN  SET  1 

SENSITIVITY  OF  BAYES  FORECAST  SCHEMES  TO  6  ASSUMPTIONS 
A  =  5;  a  =  3/4;  f  =  .8 
Average  Discounted  Cost  Per  Period 
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TABLE  3 


RUN  SET  1 

COMPARISONS  OF  PERCENT  DIFFERENCES  WITHIN  A  SCHEME 
A  =  5 ;  a  =  3/4;  f  =  .8 


TABLE  4 


RUN  SET  1 


PAIRED  t  STATISTICS  FOR  DIFFERENCE  IN  6 
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*  Significant  at  5%  level 
**  Significant  at  1%  level 
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TABLE  12 


RUN  SET  5 
SHOCKED  DEMAND 
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TABLE  14 


**  Significant  at  1%  level 


RUN  SET  5 

COMPARISON  OF  PERCENT  DIFFERENCES  IN  AVERAGE  PERIOD  COSTS 


TABLE  16 


TABLE  18 
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rats)  under  each  achame.)  As  mentioned  above,  a  comparison  was  made 
between  the  percent  of  difference  between  scheme  pairs  in  Table  6. 

The  figures  <'.n  Table  6  were  obtained  for  each  scheme  pair  by  sub¬ 
tracting  the  average  period  cost  for  the  second  scheme  from  the 
average  period  cost  for  the  first  scheme,  and  dividing  the  result 
by  the  smaller  of  the  two  average  period  costs.  In  comparing  similar 
cases  (e.g. ,  X  ■  1  ;  a  ■  1  and  A  -  1  ;  a  ■  3/4)  a  larger  percent 
difference  generally  results  with  u  »  1  than  with  a  -  3/4  .  This 
result  does  not  occur  in  all  comparisons  between  similar  cases,  but 
whan  it  does  not  the  difference  is  quite  small.  It  iu  noted  that 
when  the  case  of  trend  is  considered  (A  «  1+.5),  large  percentage 
differences  exist  in  scheme  pairs  AB;EST,  ML;EST,  and  SB;EST.  This 
is  consistent  with  the  results  of  run  set  4. 

Tables  7,  8  and  9  give  the  results  of  run  set  3  for  stationary 

demand,  where  a  ■  1  ,  and  cost  ratios  of  .8,  .3  and  .3  were  used, 

when  A  was  varied  from  .5,  1,  5,  10  and  15.  Table  7  shows  mean 

average  costs  and  standard  deviations  for  20  replications  of  40  peri¬ 

ods  each.  Table  8  displays  the  results  of  the  paired  t  tests  and 
Table  9  gives  the  percent  difference  for  those  cases  in  which  signif¬ 
icance  was  indicated  by  the  paired  t  test. 

No  significant  difference  was  indicated  among  the  ML,  SB,  and 
AB  schemes,  schemes  for  which  the  theoretical  basis  is  predicated  on 
stationary.  These  schemes  performed  significantly  better  chan  ES  and 
EST,  although  with  the  exceptio  .  of  the  low  mean  demand  (A  *  0.5  and 
A  ■  1)  cases,  the  average  percent  differences  between  ES  and  the  "sta¬ 
tionary"  schemes  were  under  5%,  It  also  appears  that  the  critical 
cost  ratio,  f  ,  has  no  particular  effect  on  the  percent  differences. 

To  check  out  a  high  value  of  f  which  in  practice  may  be  more  realis¬ 
tic,  the  cases  where  f  «=  0.95,  A  ■  0.5,  1  and  10  respectively,  were 
also  investigated.  The  results  are  presented  in  Table  16  and  are  in 
line  with  the  above  statements,  generally,  with  percent  differences 
similar  to  the  f  *  0.8  cases.  As  expected,  EST  performed  signifi¬ 
cantly  poorer  than  ES  on  the  stationary  data;  however,  the  percent 
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differences  are  generally  under  52  except  when  \  m  .5.  Even  there, 
except  for  the  very  small  f  value,  the  percent  differences  were 
under  7%. 

Thus,  when  considering  long  run  average  costs  per  period 
when  some  prior  data  exists,  if  the  demand  is  stationary,  one  is 
better  off  using  ML  or  the  Bayes  scheme  especially  for  low  demand, 
although  using  exponential  smoothing  for  moderate  or  high  demand  is 
not  too  disadvantageous.  (It  should  be  noted  that  changing  the  smooth¬ 
ing  constant (s)  used  in  ES  and  EST  might  improve  tueir  performance. 
This  effect  was  not  studied  here.) 

The  next  situations  investigated  were  those  with  a  constant 
increasing  long  term  trend.  These  results  are  presented  in  Tables 
10,  11,  and  12.  Trends  of  0.05,  0.1  and  0.5  per  period  were 
studied,  starting  with  a  base  1  of  1.  Although  the  t  test  indi¬ 
cated  significance  almost  everywhere,  when  looking  at  the  average 
percent  differences,  sizable  values  were  found  only  between  the 
"stationary"  schemes  (ML,  SB,  aB)  and  the  exponential  smoothing 
routines,  with  the  exponential  smooching  routines  showing  signifi¬ 
cantly  better  performance.  (The  maximum  percent  differences  among 
ML,  SB  and  AB  were  less  than  3%.)  For  the  trend  of  0.5  per  period, 
percent  differences  between  ES  and  the  ML,  SB,  AB  schemes  were  in 
the  70' s  when  f  **  0.8  .  A  run  in  which  f  was  set  at  0.95  shows 
percent  differences  in  the  200's  (see  Table  16).  It  appears  the 
smallest  differences  occur  at  f  =  0.5  ,  with  differences  becoming 
larger  as  f  is  increased  or  decreased,  particularly  for  f  increas¬ 
ing.  One  rather  unexpected  development  did  occur  in  this  run  set. 

In  almost  all  cases  EST  did  worse  than  ES  which  is  certainly 
counter-intuitive.  Only  when  \  =  1  +  0.5  for  f  =  0.8  was  EST 
better.  Although  the  t  test  showed  significant  differences,  the 
average,  percent  differences  were  small  (always  under  52).  It  is 
believed  the  reason  for  this  can  be  explained  as  follows. 

When  using  the  Poisson  distributions,  unless  the 

trend  is  quite  sizable,  the  period  to  period  variation  is  mistakenly 
picked  up  as  trend  in  the  EST  scheme  and  thereby  does  more  harm  than 
good  in  adjusting  the  forecast.  Figure  1  shows  a  plot  of  typically 
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generated  data  from  a  Poisson  where  the  trend  per  period  is  0.1  and 
1.0  respectively,  illustrating  the  period  to  period  variation.  For 
the  10  periods  plotted,  it  is  extremely  difficult  to  spot  trend  when 
the  trend  is  0.1  per  period. 

In  summing  up  for  the  trend  cases,  one  is  definitely  better 
off  using  exponential  smoothing  rather  than  maximum  likelihood  or 
Bayes  if  a  trend  exists.  Further,  if  the  trend  is  sizable,  EST  is 
preferred;  however,  if  the  trend  is  small,  EST  performs  slightly 
worse  than  ES. 

The  next  set  of  runs  investigated  "shock"  changes  in  demand. 

For  f  ■  0.8,  0.5  and  0.3,  A  was  changed  from  1  to  5,  1  to  10, 

10  to  1,  and  1  to  5  to  1  during  the  course  of  the  runs  (40  periods). 
Tables  13,  14  and  15  present  the  results.  For  each  of  the  shock 
runs,  the  teat  periods  are  divided  into  two  equal  segments  (except 
the  last  which  is  divided  into  3)  with  the  changes  in  A  taking 
place  between  segments. 

Once  again,  although  some  t  tests  showed  significance  among 
ML,  AB  and  SB,  the  greatest  percent  difference  was  found  to  be  0.45%; 
hence,  it  appears  that  there  are  no  "practical"  differences  among 
these  schemes.  Thus,  only  percent  differences  between  these  schemes 
and  the  exponential  smoothing  schemes  are  given  in  Table  15.  There 
appears  to  be  a  definite  advantage  in  using  ES  or  EST  over  AB,  SB 
or  ML  when  shock  changes  take  place.  Once  again  ES  outperformed  EST 
but  all  percent  differences  were  under  5%.  Results  for  f  ■  0.95 
are  shown  in  Table  16  and  show  that  the  differences  are  even  greater 
for  large  f  .  Further,  the  f  “  0.95  case  shown  in  Table  16  along 
with  those  of  Table  15  indicate  that  for  large  f  values,  the  percent 
differences  are  more  acute  when  the  mean  is  shocked  upward  while  for 
low  f  values  the  opposite  is  true.  For  the  double  shocked  case, 
differences  appear  to  be  more  pronounced  for  high  f  ;  however,  had 
the  double  shock  been  5  to  1  to  5,  this  situation  may  have  been 
reversed . 
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Th«  naxt  aerlas  of  runa  vara  concarnad  with  comparing  ML,  SB 
and  AB  for  "amall"  mean  demand  under  stationary  conditions.  Table  17 

shows  long  range  average  costs  for  X  -  0.5  and  0  -  0.05,  0.5  and 

5.0.  Once  again,  as  was  shown  in  Table  2,  the  average  cost  rate  is 

not  very  sensitive  to  0  ,  nor  to  which  of  the  "stationary"  schemes 

used  (although  if  9  greatly  underestimates  X  ,  ML  appears  somewhat 
better).  It  does  appear,  however,  it  is  preferable  to  overestimate 
Che  prior  mean  rather  than  underestimate  it,  when  f  is  "high." 

Up  to  now,  all  runs  made  had  a  pre-initialization  period  of 
5  periods  and  an  initialization  period  of  10  periods  pi  lor  to  con¬ 
sidering  costs.  To  Investigate  what  happens  in  the  Initial  periods 
without  prior  data,  a  set  of  runs  was  performed  with  no  pre- 
initialization  or  initialization,  run  length  of  first  5  and  then  10 
periods,  40  replications  and  X  -  0.5  for  6  -  0.05,  0.5  and  5.0 
respectively  at  f  *  0.8  .  It  is  summarized  in  Table  18.  We  see 
that  axcapt  where  8  is  much  greater  than  X  (6  ■  5)  there  is  essen¬ 
tially  no  difference  between  AB  and  SB.  When  6  ■  5  ,  AB  is  the 
poorer.  This  is  due  to  AB  yielding  a  greater  S*  in  the  first  peri¬ 
od  (S*(AB)  -  8,  S*(SB)  -  7  ;  see  TaLle  19).  Since  in  most  of  the 
replications,  the  total  demand  over  the  entire  10  periods  was  less 
than  7,  no  stock  outs  ever  occurred  and  AB  always  had  a  higher 
holding  cost.  Thus,  it  would  appear  it  is  not  of  any  benefit  to  use 
AB  over  SB  and  might  be  a  detriment. 

Except  for  the  0  ■  5  cases,  however,  the  Bayes  schemes  seem 
to  offer  some  advantage  over  ML  when  considering  the  first  5  periods. 
Significance  was  indicated  only  for  the  6-0.5  case;  however,  for 
10  periods,  once  again  the  Bayes  schemes  were  significantly  better 
for  0-0.5  ,  but  the  percent  differences  were  less  than  those  for 
5  periods.  For  0-0.05  ,  ML  showed  better,  although  results  were 
not  highly  significant.  Long  run  results  (Run  Set  3,  Tables  7,  8 
and  9  and  Run  Set  7,  Table  17)  showed  no  differences  among  ML,  AB 
and  SB  for  6  -  0.5  ,  but  ML  superior  for  0  -  0.05  .  Thus,  it 
appears  that  in  the  very  short  term,  Bayes  schemes  may  offer  an  ad¬ 
vantage,  especially  if  the  prior  mean  is  close  to  the  true  mean. 


28  - 


T-261 


In  order  to  batter  understand  the  similarities  and  differences 
among  ML,  SB  and  AB,  formulas  leading  to  the  stockage  levels,  S*  ,  for 
each  scheme  are  presented  below. 

t  +6 


S* 

ML:  I  Cl/x!)  e 

x-0 


S* 

SB:  l  (1/xi)  e  n+1/0 

x-0 


L  TU  g  . 

-S_-  A  +  e\x 
n+1  f_n - )  .  f 

\” +  y 

tn+*~  A  +  1  \x 
n+l/e  [J _ ,]  . 

\n  +  l/Qj 


Jo  i)  (»  +  1/8  +  1)  " 


Comparing  ML  and  SB,  we  have  already  mentioned  that  as  n  gets 

large,  then  both  converge  to 

S*  , 

l  (1/xl )  e"Cn/n(t_/n)x  -  f  . 
x**0 

Further,  for  n  ■  0  ,  they  both  yield  identical  S*  valueB.  Also, 
if  0  *  1  ,  they  are  identical.  Their  greatest  difference  would 
occur,  then,  for  the  early  periods  after  the  initial  forecast  and 
where  6  is  different  from  1.  Thus,  in  most  cases,  one  would  ex¬ 
pect  rather  similar  results  from  ML  and  SB,  especially  if  results 
are  considered  over  a  long  planning  horizon. 


In  comparing  SB  and  Ao,  it  can  be  shown  that  for  large  n  , 
these  two  schemes  should  produce  similar  results.  Letting  m  ■> 
n  +  1/6  +  1  and  r  ■  tn  +  1  ,  we  can  write 


ABi  £  (T1)*1 


s* 

Z  g(x)  -  f 
x-0 


The  negative  binomial  g(x)  has  a  moment  generating  function  (MGF) 
of 

1\  r 

MGF(t) 
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Lotting  y  ■  r/m  ,  we  have 


MGF(t) 


Lattipg  r  -*■  “  yields 


HGF(t) 


e-p/e-Pet  _  eP(et-l) 


which  is  the  MGF  for  a  Poisson,  mean  u  -  When  r  is  large,  m  must 
also  be  large  which  results  only  when  n  is  large.  Hence,  for  n 
large,  the  negative  binomial  goes  to  a  roisson  with  mean 

t  +  1 

"  ■  r/B  ■  r+i7e~  ' 

Further,  if  n  is  large,  then 

t  +  1 

p  *  r+T7? 


which  Is  the  mean  for  the  Poisson  when  using  scheme  SB.  Thus,  when 
n  is  large,  one  would  not  expect  Ve.*y  much  difference  among  ML,  SB 
and  AB.  For  small  n  ,  AB  will  tend  to  produce  higher  S*  values, 
especially  if  0  is  large.  If  0  is  small,  even  for  small  n  ,  AB 
and  SB  will  tend  to  produce  similar  S*'s  .  Table  19  compares  the 
CDF's  for  n  ■  0  ,  6  ■  0.05,  0.5  and  5.  Only  for  8  ■  5,  do  different 
S*'s  result  when  f  -  0.8  . 


VI.  CONCLUSIONS 

Keeping  In  mind  the  underlying  assumptions  and  range  of  this 
study  (namely  that  demand  per  period  is,  Poisson  with  unknown  mean  A 
and  the  cases  studied  generally  cover  A's  from  0.5  to  15),  we 
summarize  the  conclusions  as  follows. 

1.  When  considering  long  range  average  costs  per  period  with 
some  amount  of  prior  data  available,  if  the  demand  process  is  not 
stationary  (either  trend  or  shock  change  is  present)  exponential 
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smoothing  shows  suparlor  performance  over  maximum  likelihood  or 
Bayes.  If  a  sizable  linear  trend  exists,  EST  shows  up  somewhat 
better  than  ES.  Otherwise,  ES  is  preferred. 

2.  If  the  demand  process  is  stationary,  the  exponential 
smoothing  schemes  perform  poorer  chan  maximum  likelihood  and  Bayes, 
especially  for  the  lower  values  of  mean  demand  (A  •  0.5  and  1). 

For  higher  A's  however,  the  differences  were  less  than  5%. 

3.  If  no  knowledge  is  known  about  whether  the  demand  process 
is  likely  to  be  stationary,  exponential  smoothing  appears  to  be  a 
good  hedge  against  this  type  of  uncertainty. 

4.  In  no  cases  did  there  appear  to  be  substantial  differ¬ 
ences  among  ML,  SB  and  AB.  Since  ML  is  the  simplest  and  most 
understood,  it  is  recommended  in  the  case  of  stationary  demand  with 
low  mean. 


5.  When  no  prior  data  exists,  and  the  initial  period  cost9 
are  considered  for  Che  case  of  stationary  demand,  SB  appears  to  offer 
some  advantage  for  the  initial  periods,  the  advantage  being  greater 
when  6  is  close  to  A  .  Thus  a  possible  strategy  might  be  to  use 
SB  for  the  initial  few  periods  then  switch  to  ML. 

Two  further  avenues  of  exploration  would  be  of  interest. 

First,  a  continuation  of  the  initial  period  investigation  for  differ¬ 
ences  among  ML,  SB  and  AB  for  small  mean  stationary  demands  (even 
smaller  than  those  studied  here).  It  has  been  claimed  that  under 
these  conditions,  Bayes  schemes  should  show  superior  performance. 
Although  the  results  of  Table  18  are  to  some  degree  inconclusive, 
they  tend  to  indicate  this  might  be  so.  Initial  period  studies  of 
all  schemes  for  nonstationary  conditions  would  also  be  of  interest. 

Second,  it  would  be  enlightening  to  investigate  the  robust¬ 
ness  of  forecasting  schemes  if  demand  were  not  really  Poisson.  Schemes 
such  as  ES  and  ML  do  not  have  to  be  tied  to  any  assumption  concerning 


-  31  - 


T-261 


distribution  of  demand.  Often  in  practice,  what  is  done  is  to  merely 
assume  that  forecast  errors  are  normally  distributed  (this  may' or 

I 

may  not  be  true) ,  keep  a  running  calculation  of  the  sample  standard 
deviation  of  the  forecast  errors  and  use,  then,  in  Equation  (l)a 
normal  CDF  with  mean  equal  to  the  forecaat  value  and  standard  devi¬ 
ation  equal  to  the  sample  standard  deviation  of  forecast  errors. 

It  would  be  interesting  to  investigate  this  procedure  under  varying 
conditions . 

,  i 
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